function gradient = gaussCovarianceGradient(mu, sigma, x)

inv_sigma = inv(sigma);
gradient = 0.5*gauss(mu, sigma, x)*(inv_sigma*(x-mu)'*(x-mu)*inv_sigma-inv_sigma);
